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ABSTRACT 

Here we report on observations of the radio magnetar PSR J 1622^-950 at frequencies from 
1.4 to 17 GHz. We show that although its flux density is varying up to a factor of ~10 within a 
few days, it has on average decreased by a factor of 2 over the last 700 days. At the same time, 
timing analysis indicates a trend of decreasing spin-down rate over our entire data set, again 
of about a factor of 2 over 700 days, but also an erratic variability in the spin-down rate within 
this time span. Integrated pulse profiles are often close to 100 per cent linearly polarized, but 
large variations in both the profile shape and fractional polarization are regularly observed. 
Furthermore, the behaviour of the position angle of the linear polarization is very complex 
- offsets in both the absolute position angle and the phase of the position angle sweep are 
often seen and the occasional presence of orthogonal mode jumps further complicates the 
picture. However, model fitting indicates that the magnetic and rotation axes are close to 
aligned. Finally, a single pulse analysis has been carried out at four observing frequencies, 
demonstrating that the wide pulse profile is built up of narrow spikes of emission, with widths 
that scale inversely with observing frequency. All three of the known radio magnetars seem 
to have similar characteristics, with highly polarized emission, time-variable flux density and 
pulse profiles, and with spectral indices close to zero. 
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1 INTRODUCTION 

PSR J 1622^1950 was discovered in the High Time Resolution Uni- 
verse survey for pulsars and fast transients (Keith et al. 2010) cur- 
rently underway at the Parkes and Effelsberg radio telescopes. The 
pulsar's many similarities with the two previously known magne- 
tars that emit radio pulsations have placed this source in the fast 
growing group of magnetars. Magnetars are commonly thought to 
be rotating neutron stars that in addition to their emission of pulsat- 
ing radiation also undergo large bursts and outbursts of radiation in 
the X-ray and 7-ray bands (for more detailed reviews on magnetars 
see e.g. Mereghetti (2008); Rea & Esposito (2011)). The magne- 



tar group is built up of two subgroups: anomalous X-ray pulsars 
(AXPs) and soft gamma-ray repeaters (SGRs). There is however 
no longer a strict division between the two classes, as new observa- 
tions have shown that some of the sources simultaneously exhibit 
properties originally thought to belong exclusively to only one of 
the two classes (Gavriil et al. 2002; Mereghetti et al. 2009; Rea et 
al. 2009). The term magnetar originates in the sources' extremely 
high inferred surface magnetic fields (typically > 10 14 G) and it is 
believed that their radiation is powered by the energy stored in the 
magnetic fields (Duncan & Thompson 1992) instead of by the spin- 
down as is the case for ordinary pulsars. Recently, a new magne- 
tar (SGR 04 18+5729) with a considerably lower surface magnetic 
field (B < 7.5 x 10 12 G) was discovered (Rea et al. 2010), casting 
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doubts on the assumption that a high surface dipolar magnetic field 
strength is a requirement for magnetar-like activity. 

Radio pulsations from a magnetar were first detected in 2006 
from the source XTEJ1810-197 (Camilo et al. 2006), and since 
then only two other sources have been found to belong to the 
group of radio-emitting magnetars: IE 1547.0-5408 (Camilo et al. 
2007c) and PSR J1622-4950 (Levin et al. 2010). XTEJ1810-197 
and IE 1547.0-5408 are both so-called transient magnetars, that 
occasionally undergo large outbursts of X-ray emission. The ra- 
dio properties of these two pulsars have been described in detail 
in a series of papers (e.g. Camilo et al. 2007a, 2008; Kramer et al. 
2007; Lazaridis et al. 2008), which reported on features that make 
the radio magnetars stand out from the ordinary pulsar population. 
In addition to long pulse periods and high surface magnetic field 
strengths, these features include highly variable radio flux densi- 
ties, changing pulse profiles on short time scales, large amounts of 
timing noise and a flat radio spectrum. 

Both sources emit nearly 100% linearly polarized radiation at 
a large range of observing frequencies (Kramer et al. 2007; Camilo 
et al. 2007a, 2008). Analyses of the linear polarization position an- 
gle (PA) show a preferred neutron star geometry for IE 1547.0- 
5408 where the rotation and magnetic axes are close to aligned 
(Camilo et al. 2008), which at the time was supported by a low 
pulsed fraction in the X-ray (Halpern et al. 2008). Since then, X- 
ray monitoring has shown much higher pulse fraction values (Israel 
et al. 2010) and it has been shown that the low pulsed fraction in 
X-ray observed during high flux states may be due to a dust scatter- 
ing halo (Israel et al. 2010; Bernardini et al. 201 1 ; Scholz & Kaspi 
2011). This is in conflict with the aligned geometry derived from 
the radio emission of IE 1547.0-5408. The radio emission geom- 
etry analysis for XTEJ1 8 10-197 has proven difficult and different 
groups have reported on different results. Camilo et al. (2007a) re- 
port on two possible solutions for the geometry: Either the magnetic 
and rotation axes are nearly aligned, or the emission originates high 
above the surface of the star. Kramer et al. (2007) on the other hand 
derive a geometry where two emission cones must be present in the 
neutron star magnetic field. This is interpreted as either an offset 
dipole or a non-dipolar field configuration. Analyses of the X-ray 
data from XTEJ1810-197 seem to favor a non-aligned geometry 
for this source (Perna & Gotthelf 2008). 

Single pulse studies of XTE J1810— 197 are covered in great 
detail in Serylak et al. (2009). They show that the integrated pulse 
profile consists of strong spiky sub-pulses, with an overall high 
modulation index that varies between components of the pulse. 
Their analysis concludes that the radio emission from XTE J1810- 
197 is clearly different to that from ordinary pulsars. 

The radio magnetar PSR Jl 622-4950 was discovered by Levin 
et al. (2010). That paper reports on the high variability of this 
pulsar in the shape of the integrated pulse profile and in radio 
flux density on time scales of hours. The X-ray counterpart of the 
pulsar is identified, with an observed X-ray luminosity Lx(0.3- 
10 keV) « 2.5 x 10 33 ergss~\ which is in the lower end of 
the range of X-ray luminosities observed for quiescent magnetars 
(1.8xl0 33 < L x (l-10keV)< 1.2x 10 36 ; Rea&Esposito (2011); 
or see the McGill SGR/AXP Online Catalogue 1 ). The same pa- 
per also mentions that PSR J1622-4950 has highly linearly polar- 
ized radio emission and an inverted radio spectrum. Further studies 
of the radio spectrum of this magnetar have been carried out by 
Keith et al. (2011), who observed it at bands centered at 17 GHz 
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and 24 GHz, and concluded that the spectral index is close to zero 
when the flux densities from these observations are added in the 
radio spectrum calculation. 

In this paper, we will present the continued observations and 
analysis carried out for PSR J1622^1950 with the Parkes radio tele- 
scope. We will focus on four different aspects of the emission: first 
we will have a look at the flux density evolution over ~2 years of 
observations to see if the source is still as variable in flux density as 
has been measured previously and if there are any trends in how it is 
varying overall. The second part will treat the timing analysis of the 
source, with the complications that a highly varying pulse profile 
introduces to this process. We will investigate if the frequency (v) 
and frequency derivative (v) that we observe demonstrate the true 
spin-down of the pulsar or if their changes are artifacts of the pulse 
profile variations. Thereafter we will report on the polarimetry of 
the integrated pulse profiles and compare the results to the two pre- 
viously known radio magnetars. We will also attempt to derive the 
geometry of the radio emitting regions of the neutron star. Finally 
we will report on an analysis of single pulses from the magnetar 
at several observing frequencies. Throughout the paper we will re- 
view the similarities and differences between PSR J1622-4950 and 
IE 1547.0-5408, XTEJ1810-197 and ordinary pulsars. This dis- 
cussion will be incorporated in each pulsar property section in the 
paper, and we will finish by summarizing our findings in the last 
section. 



2 OBSERVATIONS AND ANALYSIS 

All data used for this analysis were collected with the 64-m dish at 
the Parkes Radio Telescope using different receivers and backends. 
Observations were made at frequency bands centered at: 1.4 GHz 
using the center beam of the Multibeam Receiver (Staveley-Smith 
et al. 1996), 3.1GHz using the '10-/50-cm' receiver and 17 and 
24 GHz using the '13-mm' receiver. 

The Parkes Digital Filterbank System (PDFB3) used to cre- 
ate the folded profiles first converts the analogue voltages from 
each polarization channel of the linear feeds into digital signals. 
It then produces 1024 polyphase filterbank frequency channels that 
are folded at the apparent topocentric period of the pulsar into 1024 
pulsar phase bins, and written to disk every 20 s. Four Stokes pa- 
rameters are recorded. To determine the relative gain of the two 
polarization channels and the phase between them, a calibration 
signal is injected at an angle of 45° to the feed probes. The data are 
analysed off-line using the PSRCHIVE package 2 (Hotan et al. 2004) 
and corrected for parallactic angle and the orientation of the feed. 
The position angles are also corrected for Faraday rotation through 
the interstellar medium using the nominal rotation measure. 

For the single pulse analysis the baseband data recording and 
processing system known as the ATNF Parkes Swinburne Recorder 
(APSR; van Straten & Bailes 201 1) was used. 



3 RADIO LIGHT CURVE 

As reported in Levin et al. (2010) the 1.4 GHz flux density of the 
integrated pulse profile varies greatly between observations. Since 
then we have collected about one more year of data on this source. 
Plotted in the bottom panel of Fig. 1 is the integrated flux density 
of each observation made at 1 .4 GHz at Parkes since the discovery 

2 See http://psrchive.sourceforge.net 
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Figure 1. Top plot: Variation of frequency derivative with time. The error bars in the x-range indicate the time span of the included TOAs. The large variations 
in the value of the frequency derivative indicate changes in absolute value and sometimes also in sign of the second derivative. In some cases we have included 
a fit to a shorter data set, even if all the TOAs in that set are already included in a fit of a longer data set, to show the direction of the frequency derivative 
change. In total, the values of the frequency derivative varies up to a factor of ~2, Bottom plot: Variation of flux density at 1 .4 GHz with time. The mean and 
RMS values were calculated using the same sets of observations as were used in the frequency derivative analysis. 



in April 2009. The new data points begin around MJD=55240. It is 
clear from this plot that the peak observed flux density during the 
last ~200 days is only about half the value observed during the first 
~200 days after the discovery, suggesting an intrinsic long-term 
decay of the flux density. However, it is possible that the magnetar 
is still just as variable and with as high peak flux density as before 
but that the higher flux density points are missed during the later 
time span due to the smaller sample. Simply fitting a line to the data 
points results in a slight slope, giving a decline of the average flux 
density of ~2 mJy for the 700 days of observing. To more easily 
visualise the flux density decline we have divided the data points 
into smaller sets of about 30 days each and calculated the mean 
flux density and the root mean square (RMS) for each of these data 
sets. The time span for each data set is the same as we used for 
calculating different rotation frequency derivatives as described in 
Sec. 4 below. The mean and RMS of each data set is plotted on top 
of the data points in the bottom panel of Fig. 1. 

Anderson et al. (2011b) report on recent observations of 
PSR J1622^950 with the Chandra X-ray Observatory and the Aus- 
tralia Telescope Compact Array (ATCA) that were collected within 
the framework for the "ChlcAGO" project (Anderson et al. 201 la). 
They observed the magnetar with the ATCA simultaneously at fre- 
quency bands centered at 5.0 GHz and 9.0 GHz on November 22, 
2008 and December 5, 2008. The flux densities are 33.0±0.3mJy 
and40.4±0.3 mJy at 5.0GHz and 30.9±0.6mJy and 31.9±0.6mJy 
at 9.0 GHz for the two observation sessions respectively. This is 
significantly higher flux densities than the values measured on De- 
cember 8, 2009 and February 27, 2010 with the ATCA published 



in Levin et al. (2010): 13±lmJy at 5.0GHz and 14.3±0.8mJy at 
9.0 GHz. These values indicate a decrease in flux density of about 
68% at 5.0 GHz and 55% at 9.0 GHz over one year (Anderson et 
al. 2011b). With this in mind it seems likely that the flux density 
decline observed recently at 1 .4 GHz with Parkes is a real intrinsic 
decay of flux density in the magnetar. 



4 TIMING 

To analyse the rotational history of a pulsar, it is conventional to 
integrate over a set number of rotations, enough to create a stable 
pulse profile for each observation. These profiles are then aligned 
with the help of a standard profile, unique to that particular pul- 
sar, to create a list of the pulse times of arrival (TOAs) at the tele- 
scope. The list of TOAs is used to determine a more accurate period 
and spin-down rate of the pulsar as well as its precise position in 
the sky, a more precise value of its dispersion measure (DM), etc. 
When performed on data from a large enough time span, using this 
method generally results in values for the pulsar parameters with 
very high accuracy. 

However, due to the variability of pulse profile at 1 .4 GHz, the 
timing analysis of PSR J 1622-4950 is more complex than for an 
ordinary pulsar. Here we have employed the same tools in the tim- 
ing procedure as were used in Levin et al. (2010), i.e. TOAs for the 
pulses were calculated by using a model that describes the differ- 
ent components of the profile, rather than using a standard profile. 
The model was created using software from the PSRCHIVE pack- 
age (Hotan et al. 2004) by fitting scaled von Mises functions (von 
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Figure 2. Timing model. The dashed lines show the von Mises functions 
used to build up the profile, which added together give the solid line that 
follows the pulse shape. The bottom noise level shows the emission that is 
left after subtracting the model. Profile from 29 June 2009. 

Mises 1918) to the pulse profile of one of the observations where 
all the components were present. This profile and the correspond- 
ing model is shown in Fig. 2. The TOAs are then created by letting 
the amplitudes of the components vary but keeping the separations 
fixed, while fitting the model to each observation. The timing anal- 
ysis was made using the TEMPO software 3 . 

However, we cannot obtain a coherent timing solution for our 
full data span even using this method. Instead we have looked at 
data from shorter time spans and fitted the observed frequency (v) 
and frequency derivative (i>) for each set of timing points sepa- 
rately, in an attempt to quantify how much and in which direction 
the true spin-down is changing with time. The number of points 
included in each set is dependent on how long we could get a sta- 
ble i> with a reasonable error (largest error in Fig. 1 is i> err = 0.21 
x 10~ 3 s~ ). For this timing analysis the position is held fixed 
at that constrained by the X-ray counterpart: R.A. = 16:22:44.80, 
Dec. = -49:50:54.4 (Levin et al. 2010) and the DM is set to 820 
cm _3 pc. The result of this analysis is shown in the upper panel of 
Fig 1 . It is evident from this plot that ii has been changing with a 
factor of ~2 since the discovery, as was also stated in Levin et al. 
(2010). However, the second derivative of the frequency (v) is also 
changing rapidly in magnitude and direction with time. 

A comparison of these results to the timing analysis carried 
out previously for the other two radio magnetars, IE 1547.0-5408 
(Camilo et al. 2007b) and XTEJ1810-197 (Camilo et al. 2008) 
shows clear similarities, but also differences. In all three cases 
the v is changing greatly as a function of time. In both of the 
other sources the v seems to vary more smoothly than it does for 
PSRJ1622-4950 and they both have steady trends along a fairly 
straight line (allowing for some "wobbling" on the way). How- 
ever, the sign of the estimated v is different for the two sources: 
XTE J1810-197 has a positive i> and IE 1547.0-5408 has a nega- 
tive one. The timing analyses for both the other two magnetars were 
made using data that were collected more regularly and frequently 
than our observations, which may contribute to the smoother look- 
ing frequency derivative evolution. 



3 See http://www.atnf.csiro.au/research/pulsar/tempo/ 



Camilo et al. (2007b) quantify the timing noise by looking 
at the magnitude of the cubic term of a Taylor series expansion 
of rotational phase over a time interval t, i.e. i> t 3 /6 (Arzouma- 
nian et al. 1994). By using this expression they get about 60 cycles 
over 6 months for IE 1547.0-5408 and 20 cycles over 9 months for 
XTE J1810-197. The same calculation for PSRJT 622^1950 yields 
about 250 cycles over the entire 20 months that the pulsar has 
been observed since the discovery. If we instead look at only the 
100 days when the v is steadily increasing (from MJD~55150 to 
MJD~55250) the value is about 5 cycles over 3.3 months. 

One option that could explain the unusual behavior of v in 
PSR J1622-4950, is if the source went through a glitch shortly be- 
fore the discovery observation at MJD=54939. Glitches are more 
frequent for younger pulsars than for older sources, but in general 
unusual in ordinary pulsars (Espinoza et al. 201 1). For AXPs how- 
ever, glitches have been observed in nearly all known sources (Dib 
et al. 2008) and it would therefore be feasible to observe a glitch 
also for PSR J1622-4950. Unfortunately it is very hard to constrain 
if a glitch has occurred without any data from around or before the 
time of the possible glitch, but the possibility is worth noting for 
future timing efforts of this source. 

By comparing the two panels in Fig. 1, there is an apparent 
correlation between reduced flux density and increased v, again 
especially between MJD « 55100 and MJD « 55300. Calculat- 
ing Pearson's correlation coefficient, ri x ^ y ) =cov(x,y)/a x a y (where 
cov(x,y) is the covariance of x and y, and a is the standard devia- 
tion), for correlation between v and the mean values for each time 
span in Fig. 1 gives r = -0.73, which corresponds to a 3 sigma cor- 
relation. It is however hard to constrain that these two features are 
direct consequences of each other, and we stress that care should 
be taken not to infer too much from these correlation results. 

Beloborodov (2009) describes a model in which starquakes 
or glitches in the neutron star cause the magnetic field lines to 
twist. The author describes the electrodynamics of the untwisting 
of the field lines as they relax back towards their initial state. He 
considers this model in the context of the magnetars, in particu- 
lar XTE J1810-197, where the reduction in the radio flux density 
over time, the changes to the torque and the X-ray observations all 
conform roughly with his calculations. This model may also be ap- 
plicable here as we see a similar decrease in flux density over time 
and large variations to the torque although we are hampered by a 
lack of information in the X-ray band. It is also hard to see how 
the gradual untwisting of magnetic field lines can cause the rapid 
profile and polarization variations that we see as these appear to 
oscillate back and forwards between states on a time-scale much 
shorter than the decay of the radio flux density. 



5 POLARIMETRY 

We have collected polarimetric data for PSR J1622-4950 with the 
Parkes telescope at four different observing frequencies, with bands 
centered at: 1.4 GHz (the centre beam of the Multibeam receiver), 
3.1 GHz (10-/50-cm receiver) and 17 and 24 GHz (13 mm receiver). 

The high frequency observations (at 17 and 24 GHz) were re- 
ported on in Keith et al. (2011), and we note some striking simi- 
larities in the polarized emission over all observed frequencies. As 
well as having a very wide (~50%) duty cycle for the integrated 
profiles, the linear polarization is often close to 100% of the to- 
tal intensity. However, similar to the flux density and pulse profile 
shape, also the polarization of the integrated profiles is changing 
between observations. 
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Figure 3. Examples of observations at 1 .4 GHz from the four different polarization groups described in Sec 5 . 1 . Category I - IV are shown from top left to 
bottom right. Total intensity is shown in black, linear polarization in red and circular polarization in blue. The integration time for each observation is noted in 
the panels, however there are both shorter (^180 s) and longer (~600 s) observations present in all categories. Note the short time scale of the variations: all 
four observations were collected within the same calender month (December 2009). 



5.1 Polarisation categories 

In an attempt to get an overview of the way the polarization is vary- 
ing we have collected observations with similar characteristics in 
groups, giving us four separate categories with different properties. 
Category I is the most frequent mode, with almost half (48%) of 
all categorised observations. It collects the observations where the 
linear polarization is > 50% of the total intensity, the circular po- 
larization is low and the position angle (PA) has a steep and con- 
sistent swing. In category II the linear polarization is much lower 
throughout the profile. In the cases where a second component is 
present, the fraction of linear polarization is higher in the trailing 
edge than in the leading edge of the profile. Similar to category I, 
the circular polarization is low and the PA swing is steep and con- 
sistent. 24% of our observations belong to this group. For category 
III (1 1% of the observations) the most prominent feature is the shal- 
low PA curve. In addition the profiles also tend to have low linear 
polarization in the leading edge of the profile and a higher value of 
the circular polarization than for categories I and II. The last 17% 
of the observations do not fit into any of the first three groups, and 
hence will end up in category IV. Here we have collected the obser- 
vations with jumps and other irregularities in the PA swing, change 
of handedness in the circular polarization and large changes in lin- 
ear polarization within the pulse profile. Examples from all four 
groups are shown in Fig 3. 

When looking at the time evolution of the polarized emission, 
by analysing observations from the different groups in time order, it 
does not seem like the variations are following any preferred order, 



but are fairly random in time. The large variations in linear polariza- 
tion and PA for PSR J 1622-4950 are at odds with what is seen for 
XTEJ1810-197 by Kramer et al. (2007). They observe an evolu- 
tion in PA swing over a time-scale of weeks, but very few changes 
on shorter time-scales. Also Camilo et al. (2007a) observed that 
the general polarization properties of XTEJ18 10-197 do not seem 
to vary with time as the total intensity changes. That is, the lin- 
ear polarization is always close to 100% of the total intensity, and 
the circular polarization component is low. This suggests that the 
observed profile variations are not caused by changes in the mag- 
netic field geometry of the emission regions for this source. Also 
in the case of IE 1547.0-5408 there is little variation in linear po- 
larization and in measured PA swing with observing frequency and 
time (Camilo et al. 2008). The circular polarization however in- 
creases with decreasing frequency, and is overall higher than for 
XTEJ1810-197. 



5.2 Rotating Vector Model predictions on the neutron star 
geometry 

By analysing the linear polarization and how its position angle is 
varying across the pulse profile, predictions on the angles of the 
rotation and dipole axes can be made. The rotating vector model 
(RVM; Radhakrishnan & Cooke 1969) states that the pulsar emis- 
sion beam has its base close to the dipole axis of the pulsar mag- 
netic field and is observed through rapid swings of the linear po- 
larization position angle over the pulse phase. When using the 
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Table 1. RVM fits to 3.1 GHz data. All angles are given in degrees. 



Obs ID 


MJD 


x 2 




C 


a 


00 


s091230_210826 


55195.88 


1.47 


-5.7±13.6 


15.1±28.8 


36.0±60.2 


141.6±9.7 


s091230_215432 


55195.91 


1.28 


-11.6±3.4 


15.5±9.8 


28.6±17.2 


137.5±2.0 


slOOlOl-201716 


55197.85 


1.01 


-2().6±2.8 


20.7±4.0 


46.0±7.8 


149.7±1.6 


1100615.074911 


55362.33 


1.34 


-20.5±1.8 


14.7±5.8 


28.3±10.4 


137.0±1.1 


tl00825_102356 


55433.43 


1.62 


-23.7±6.1 


16.3±22.2 


25.0±33.5 


157.4±2.2 


tll()116_211759 a 


55577.89 


2.40 


-4.8±1.5 


13.2±7.0 


22.2±11.3 


192.9±0.9 


S1104KL210120 


55661.88 


1.69 


-29.6±6.4 


7.1±21.9 


20.7±58.0 


122.6±7.3 


si 10410.212455 


55661.89 


1.10 


-26.5±3.2 


9.9±9.7 


24.7±21.8 


128.9±3.3 



Notes: 

a Observation with orthogonal PA jump 



RVM there are a number of effects that are very difficult to take 
into account and hence are often ignored. These effects include 
rotational sweepback of the magnetic field lines (Dyks & Hard- 
ing 2004), propagation effects in the pulsar magnetosphere (e.g. 
Petrova 2006), emission height differences (Dyks 2008) and mul- 
tipolar components of the magnetic field. There are a number of 
papers discussing these effects on the magnetic field geometry for 
magnetars in particular, and the possibility that higher order multi- 
poles are contributing to the magnetic field structure near the mag- 
netar surface (e.g. Thompson et al. 2002; Beloborodov 2009; Rea 
et al. 2010; Turolla et al. 2011). We will disregard these effects in 
this paper, but it is important to keep in mind that the magnetic field 
topology may well deviate from the simple dipole model. 

In the case of PSR J 1622-4950, the highly varying values of 
the PAs in the 1 .4 GHz observations make it difficult to find a con- 
sistent solution to the geometry of the neutron star emission. RVM 
fits to data from different days give different answers depending on 
the parameters of the PA curve for that particular observation. The 
scatter broadening of the single pulses at 1 .4 GHz (as described in 
section 6.2) could contribute to some of the large changes in linear 
polarization and PA swing that we observe at this observing fre- 
quency, and hence we have focused the RVM fits to the less scatter 
broadened and apparently more stable 3.1 GHz data. Tab. 1 shows 
the values from the 3.1 GHz RVM fits. Even though the best fit an- 
gles vary between observations we are able to put some constraints 
from these fits. The angle between the spin axis and the pulsar- 
observer line of sight, is small in all observations with C, < 20° 
and the angle between the spin axis and the magnetic pole, a, is al- 
ways just a few degrees higher than £ with 20° < a < 46°. The re- 
sulting angle, j3, places the line of sight between the magnetic pole 
axis and the positive rotation axis with values —25° ^ ^ —8°. 
The results imply that the pulsar has close to aligned magnetic and 
rotation axes and that the line of sight remains within the emission 
beam for large parts of the rotation, which in turn is implied by the 
wide duty cycle of the integrated profiles. 

RVM predictions for IE 1547.0-5408 are described in Camilo 
et al. (2008). They carried out polarimetric studies with several dif- 
ferent receivers at the Parkes telescope at five different observing 
frequencies ranging from 1.4 to 8.4 GHz and at the Australia Tele- 
scope Compact Array (ATCA) at frequency bands centered at 18 
and 44 GHz. They report on a slow PA sweep that has an abso- 
lute value identical at all observed frequencies, and their RVM fit 
suggests nearly aligned rotation and magnetic axes. At the time, 
this result was strengthened by a low pulsed fraction in the X-ray 
(~7%) measured as the source was in quiescence, shortly before 
the 2007 outburst (Halpern et al. 2008). It has subsequently been 
suggested that the low pulsed fraction may be due to a dust scat- 



tering halo around the magnetar (Tiengo et al. 2010; Olausen et al. 
2011). Observations of the magnetar at the time of the 2008 Oc- 
tober outburst showed much higher pulsed fraction values (~20%) 
which increased to ~50% during the following few weeks (Israel 
et al. 2010). The anti-correlation between the flux and the pulsed 
fraction that has been observed for IE 1547.0-5408 (Israel et al. 
2010; Bernardini et al. 2011; Scholz & Kaspi 2011) together with 
the highly variable pulse profiles as a function of time, may sug- 
gest that the magnetospheric geometry is variable during the initial 
phases after an outburst and may not be easily related to the geome- 
try in quiescence (Israel et al. 2010). Hence, even though the higher 
pulsed fraction measured for IE 1547.0-5408 in more recent obser- 
vations conflict with the aligned geometry derived from the radio 
data, we do not see these results as strong enough to completely 
rule out an aligned geometry. 

Analyses of XTEJ1810-197 also result in different geome- 
tries for different research groups and wavelengths. Camilo et al. 
(2007a) analysed radio polarization observations of XTEJ1810- 
197 collected with the Parkes Telescope at three observing frequen- 
cies (1.4, 3.2 and 8.4GHz). Their analysis shows a shallow swing 
of the position angle of the linear polarization, with values that 
yield two possible solutions for the geometry of the magnetar. Ei- 
ther the magnetic and rotation axes are nearly aligned, or the emis- 
sion originates high above the surface of the star. Using X-ray data, 
Perna & Gotthelf (2008) also estimated the viewing geometry of 
XTE Jl 8 10-1 97. They determine the allowed minimum and max- 
imum angles between the line of sight and the emission hot spot, 
and find that the range of the minimum value is compatible with 
very small angles (including zero) while the maximum is always 
large (> 60°). It has been shown that the peaks of the radio and 
X-ray pulses from XTE Jl 8 10-197 are well matched (Camilo et al. 
2007b), which suggests that the radio emission axis and the hot spot 
axis are very close to aligned. Perna & Gotthelf (2008) found that 
the high emission height solution from Camilo et al. (2007a) was 
well compatible with their results and that an aligned geometry was 
unlikely. However, they were not able to make a formal statistical 
comparison with the results from Camilo et al. (2007a). Concur- 
rently, Kramer et al. (2007) also published results from radio po- 
larization observations of the same magnetar. Their simultaneous 
multifrequency observations (at 1.4, 4.9 and 8.4 GHz) were done 
with three European telescopes: the 76-m Lovell radio telescope in 
the UK, the 94-m equivalent Westerbork Synthesis Telescope in the 
Netherlands, and the 100-m radio telescope at Effelsberg in Ger- 
many. The main difference in their results compared to Camilo et 
al. (2007a) is that Kramer et al. (2007) include studies of the single 
pulse emission and identify an inter-pulse in addition to the main 
pulse during some of their observations and at some frequencies, 
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that was not visible in the Parkes data. The PA values in the inter- 
pulse are observed to vary with time, which complicates the geom- 
etry analysis. Instead of fitting a RVM to the entire pulse profile 
the authors analysed the main pulse and the inter-pulse separately, 
with a geometry where two emission cones must be present in the 
neutron star magnetic field as result. This is interpreted as either an 
offset dipole or a non-dipolar field configuration. 

The many different geometry analyses that have been carried 
out for the three radio magnetars make a comparison between them 
complicated. It is still a bit unclear but we can not completely rule 
out that the three sources are all aligned, as hinted by the RVM. 
We find this solution very tempting, as it could prove insight in the 
radio behavior of the magnetars. Since the probability of observing 
a pulsar is lower for a smaller a, aligned geometries of the radio 
magnetars could explain why only three of the 23 currently known 
magnetars and magnetar candidates have observed radio pulsations. 
So far no X-ray pulsations have been observed from PSR J 1622- 
4950. Detection of pulsations in the X-ray might provide additional 
constraints on the geometry of PSR J1622-4950 and would help us 
to investigate this question further. 

5.3 Pulse profile alignment 

At 3.1 GHz the amount of linear polarization is in general very 
large. All of our observations at this frequency band belong to cat- 
egory I or category IV in the classification above. Even though the 
amount of linear polarization seems fairly constant, the integrated 
pulse profile still has a varying shape. The problem with aligning 
profiles of different pulse shapes was dealt with in the timing case 
by making a model of von Mises functions (see Section 4 and Fig 
2). Seeing that the PA is often similar for all 3.1 GHz observations, 
we have also aligned the profiles by using the value of cj>o that was 
estimated from the RVM model. The result is shown in Fig 4. Com- 
paring the two alignment methods, it again becomes obvious that 
the polarization from the magnetar is changing with time. At first 
glance it might seem like the aligning after the total intensity pro- 
file (the bottom plot of Fig. 4) is the one that should be preferred. 
However, a closer look at the PAs shows not only a shift in absolute 
value of the PAs but also that the PA curve for some of the ob- 
servations have different slopes. For example, this can be seen by 
comparing the brown crosses with the dark blue ones in the upper 
panel of the bottom plot in Fig 4. Both PAs have similar slope in 
the leading edge of the pulse, but towards the trailing end the brown 
PA curve flattens out much faster than the dark blue one. The same 
slope difference is of course also present when we align the profiles 
after the value of cj>o, but by looking at the top plot of Fig 4 (and 
disregarding from the light blue curve, which is offset from the oth- 
ers) even with the different slopes we find that the PAs align fairly 
well for the different observations. The profiles however do not. 

If the rotation measure (RM) of the interstellar medium would 
change with time between observations, that could explain the off- 
set in absolute value of the PAs in the bottom panel of Fig. 4. How- 
ever, by estimating the RM value for two of the observations where 
the offset is large (e.g. for the dark blue and the green curves) we 
do not see a large enough change between our observations. 

It is hard to say which (if any) of the two alignment meth- 
ods describes the true magnetar emission. One way to get a better 
handle on this would be to use the polarization information in the 
timing of the pulsar by calculating TOAs using the <po values from 
each observation. The errors from this timing analysis could then 
be compared to the errors from the timing model described in Sec. 
4. Such an analysis would only be possible to do with 3.1 GHz (or 



higher frequency) data, since part of the linear polarization emis- 
sion at 1.4 GHz is affected by interstellar scattering (see Section 
5.4), which causes the degree of linear polarization for many of 
the 1.4 GHz observations to be too small or the PA swing to be 
flattened, preventing a reliable RVM fit. The integrated magnetar 
emission also seems to be more stable at 3.1 GHz than it does at 
1.4 GHz, which would help in the timing for both timing methods. 
Unfortunately the number of 3.1 GHz observations carried out at 
this point is not large enough and the observations that exist are not 
sampled densely enough to allow for such an analysis at the present 
time. 



5.4 Depolarization 

In general, we observe a lower degree of linear polarization at 
1.4 GHz than at the higher observing frequencies. A similar trend 
is seen in IE 1547.0-5408 by Camilo et al. (2008). This magne- 
tar is at a similar DM as PSRI1622-4950 and they are both posi- 
tioned close to the Galactic disk, which implies that both pulsars 
have a fairly high scattering timescale as predicted by the NE2001 
model (Cordes & Lazio 2002): IE 1547.0-5408 has DM = 830 ± 
50 cm~ 3 pc and (l,b) = (327.23, -0. 13) (Camilo et al. 2007c) which 
gives a scattering time scale of ~70 ms at 1 GHz compared to DM 
= 820 ± 30 cm- 3 pc at (l,b) = (333.85, -0. 10) for PSRI1622^1950 
(Levin et al. 2010) which results in ~95 ms at the same frequency. 
Camilo et al. (2008) explain this depolarization in IE 1547.0-5408 
at lower observing frequencies partly as an effect of the interstellar 
scattering of the pulse profile at these frequencies. This causes the 
PA to rotate through the different phases of the pulse, which will 
be mixed at the observer and will thus reduce the apparent linear 
polarization. They also give deviations in rotation measure (RM) 
over different paths taken by the scattered rays as a further possi- 
ble reduction effect. They conclude that scattering effects can only 
be responsible for part of the depolarization. Hence, even though 
some of the depolarization at lower frequencies for PSR J 1622- 
4950 could be intrinsic to the source, it is likely that scattering ef- 
fects are also responsible for some depolarization in our case. In 
addition, we do observe a large amount of scattering in the single 
pulses at 1.4 GHz (see Sec. 6.1), which further justifies this hypoth- 
esis. Unfortunately, due to issues with the observing system at the 
time of collection of the 1 .4 GHz single pulse data, we are not able 
to polarization calibrate this data, and hence will not be able to 
analyse which effect the depolarization has on the 1 .4 GHz single 
pulses. The single pulse data collected at an observing frequency 
of 3.1 GHz is almost 100% linearly polarized, which is discussed 
in Sec. 6.2 below. 



6 SINGLE PULSES 

Single pulse analyses have been carried out at three observing fre- 
quencies, with bands centered at 1.4, 3.1 and 17 GHz collected us- 
ing the APSR backend at Parkes. Due to the large scattering effects 
at 1.4 GHz (discussed in Sec. 6.1 below) and the poor time reso- 
lution obtained at 17 GHz (1024 bins over the pulse profile), most 
of the single pulse analysis carried out for this paper has focused 
on 3.1 GHz data. In addition, one archival observation, collected 
within the frame work of the Methanol Multibeam Survey (Bates 
et al. 2010) at Parkes using the Analogue Filterbank at a frequency 
band centered at 6.6 GHz has been analysed. 
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Figure 4. Pulse profiles and PAs for 3.1 GHz observations without orthogonal PA jumps. The flux density in the lower panels of the two plots are normalised 
after the maximum flux density of the brightest observation. The profiles of the same color refer to the same observation for both alignments. Top: The PAs 
are aligned after the value of 4>o from the fit to the RVM model. Hence the maximum value of 5ij)/54> for each observation is centered at phase 0.3 Bottom: 
The profiles are aligned after the total intensity by centering the leading edge of each profile at phase 0.3. 
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6.1 Scattering 

Due to the high DM of PSR J 1 622^950, the single pulses are likely 
broadened enough by interstellar scattering at 1 .4 GHz that it will 
affect the pulse width at our time resolution. The pulse broadening 
due to scattering at the pulsar position and distance is ~17ms at 
1.4 GHz according to the NE2001 model (Cordes & Lazio 2002), 
assuming Kolmogorov scalings from 1 .0 GHz, but the large uncer- 
tainties in the model indicate that this value could deviate by up 
to a factor of ten (as has been shown by e.g. Bhat et al. 2004). 
To better estimate the total scattering we analysed the widths of 
bright single pulses in the 3.1 GHz data by dividing the observed 
frequency band up in parts and calculating the pulse broadening 
over the bandwidth for that observation. This resulted in ~8.7 ms 
smearing over a 768 MHz band centered at 3.1 GHz. Again by as- 
suming Kolmogorov scaling (T scat ter oc v~ a , using scaling index 
a « 4.0 as an estimation for high DM pulsars (Lohmer et al. 2001; 
Bhat et al. 2004)) we calculate a scattering of ~200 ms at a band 
centered at 1 .4 GHz, which is of the same order of magnitude as 
the measured value of the single pulse widths at that observing fre- 
quency (see Sec. 6.3 and Tab. 2). This implies that single pulses at 
1 .4 GHz will be highly dominated by scattering effects and hence 
we will not include the 1.4 GHz data in the single pulse analysis. 
The same analysis results in 0.41 ms of smearing at 6.6 GHz and 
9.3 ^s at 17 GHz. 



6.2 Single pulse characteristics 

Although the integrated pulse profile for PSR J 1622-4950 is often 
very wide (with a ~50% duty cycle), when analysing each rota- 
tion of the pulsar separately, it is clear that the profiles are built up 
by the collection of much narrower pulses (see Fig. 5 and Tab 2). 
Each pulse consists of one or a few narrow spikes, at all observed 
frequency bands. 

The top right panel of Fig. 5 shows one of the brighter single 
pulses observed at 3. 1 GHz and the bottom right panel is a zoom-in 
on the two brightest components of that particular rotation. From 
this plot we can see how each spike is almost 100% linearly polar- 
ized with no or very little circular polarization. The position angle 
of linear polarization of the separate components are well resolved 
and seem to be fairly fiat across the spike. Looking at the combined 
PAs for the rotation in the top right panel, it is clear that they fol- 
low the PA swing of the total integrated profile. When we look at 
the single spike PAs in more detail, we find they are often in agree- 
ment with the total PA swing, but there are also occasions when the 
single PAs are much steeper than the integrated PA. This results in 
that we sometimes see small 'wiggles' in the PA swing that often 
correspond to emission peaks in the total intensity of the integrated 
profile. Similar wiggles in the PA swing are also seen in some ob- 
servations of XTE J1810-197 (Kramer et al. 2007). 

At 3.1GHz a histogram of the phases of the single pulses 
matches the total intensity profile well (see Fig 6). When only the 
brightest 10% of the spikes are taken into account, we get the blue 
distribution in Fig. 6. These spikes seem to be spread out over al- 
most the entire pulse profile, but with a preference for certain phase 
bins. Each rotation consists of up to 15 spikes of emission, with 2 
or 3 spikes per rotation being most frequent. The separation be- 
tween emission spikes varies up to the full integrated pulse width, 
but with a preferred value of ~ 170 ms. 



Table 2. Widths of single pulse spikes at different observing frequencies. 
The upper limit given for the 17.0 GHz pulse width is the value of two 
phase bins. The scatter broadening is measured at 3.1 GHz and scaled using 
Kolmogorov scaling, with index a = 4.0, to estimate the values at the other 
frequency bands. 



Observing Average width Scatter broadening Width of 

frequency of single pulses of single pulses integrated profiles 
[GHz] [ms] [ms] [ms] 



1.4 215 200 2163 

3.1 26 8.4 1946 

6.6 13 0.41 1298 

17.0 SJ 8.44 0.0093 1730 



6.3 Width of single emission spikes 

The width of the total integrated pulse profiles scales with ob- 
serving frequency in ordinary pulsars. This phenomenon is usually 
thought to be due to radius-to-frequency mapping, which suggests 
that the emission at different frequencies originates from differ- 
ent altitudes above the polar cap (Cordes 1978). A study of the 
integrated profile width is not as useful for PSR J 1622-4950 un- 
less simultaneous multifrequency observations are considered, as 
the widths of the profiles are changing with time when the profiles 
vary. However, in a similar way to the integrated profiles, we ex- 
pect to see a decrease in the width of single emission spikes if the 
emission at different frequencies are emitted at different altitudes. 
Indeed, there have been indications that the pulse width is scal- 
ing down with increasing observing frequency in previous pulsar 
work (see e.g. Kramer et al. (2002) and references therein). This is 
even more clearly demonstrated in PSR J 1622-4950, see Tab. 2. By 
comparing the average widths of the single pulse spikes at the three 
highest observed frequencies (3.1, 6.6 and 17 GHz) and correcting 
for the broadening due to scattering, we calculate a frequency de- 
pendence for the width of the single emission spikes for this source 
as r oc ^-°- 62 ± 012 . Kramer et al. (2002) carried out a similar 
analysis for the width of the micro-structure in PSRB1 133+16 and 
found a much flatter dependence: r M oc j,- 00 ^ 010 . 

Kramer et al. (2002) also confirmed a relationship between 
micro-structure width and pulse period first suggested by Cordes 
(1979) as a linear dependence. The single emission spike widths for 
PSR J1622-4950 at similar observing frequencies do not fit this re- 
lation, but the width at 17 GHz (^8.44 ms, an upper limit estimated 
to two time bins) is closer to the predicted value of 3.0 ms. This may 
suggest that the scattering at the lower observing frequencies is un- 
derestimated, and that the 17 GHz width is the true width of the 
single emission spikes. However, since the scattering at 3.1 GHz 
is the measured value, we find it more likely that the width scal- 
ing with observing frequency is intrinsic, which indicates that per- 
haps separate relations are required for different frequency bands. 
Other reasons to why the emission from PSR J1622-4950 does not 
fit the relation could be that the pulse structure that we observe 
in PSR J1622^1950 is not of the same origin as microstructure, or 
that magnetars follow a different relation than ordinary pulsars. The 
single pulse widths for the other two radio-emitting magnetars are 
in the same order of magnitude as for PSR J 1622-4950, and with 
the discovery of more sources it will be possible to investigate if a 
similar dependence for magnetars is valid. 
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Figure 5. Left: 3.1 GHz single pulse stack with the corresponding integrated profile on top. Linear polarization is shown in red, circular polarization in blue 
and the total intensity in black. Top right: Pulse profile of one of the brighter single pulses showing the very narrow spikes of emission. Bottom right: Zoom in 
on the two brightest spikes of emission in the rotation above. 
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Figure 6. 3.1 GHz pulse phase histogram. The red boxes show all pulses, 
while the blue boxes only account for the brightest 10% of the spikes. The 
black line is the total intensity profile averaged to match the top value of the 
histogram. 



6.4 Pulse-flux distribution 

By looking at the flux density of single spikes of emission we have 
calculated a pulse flux distribution for the single pulse spikes at 
3.1 GHz. A conventional energy distribution calculation often only 
records the peak flux density in the on-pulse region for each ro- 
tation and compares the resulting histogram to a histogram of the 
peak flux density in the off-pulse region. To account for the fre- 
quent occurrence of multiple emission spikes in a single rotation of 
PSR J1622-4950 we have approached this analysis differently. The 
on-pulse phases were chosen as all spikes that consisted of at least 
3 consecutive phase bins with a signal stronger than 3 sigma, and 
the corresponding flux value for each spike was added to the his- 
togram. The histogram in Fig. 7 was created after adding the phase 
bins to a total of 5 12 bins over the profile. The left panel shows the 
flux distribution and the right panel shows the 10 base logarithm of 
the flux. From this it is clear that the flux distribution from the mag- 
netar is closer to a log-normal rather than Gaussian distributed. It 
is also evident from these plots that most spikes are of similar flux 
density and that there were no giant pulses observed from the mag- 
netar. These properties are similar to the general ordinary pulsar 
population, for which a large fraction of the measured flux density 
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Figure 7. Pulse-flux distribution at 3.1 GHz. Left: Histogram of flux/mean 
flux of the single pulse spikes. Right: Histogram of the logarithm for the 
same data. The solid line shows the best fit of a Gaussian to the histogram, 
displaying the lognormal distribution of the flux density. 



distributions are close to log-normal and most do not emit giant 
pulses (e.g. Cairns et al. 2004; Burke-Spolaor et al. submitted). 
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Figure 8. The bottom plot shows the modulation index (dashed line) and 
the R-parameter (dotted line). For reference, the integrated pulse profile is 
plotted in the top panel. 



6.5 Pulse modulation 

To get an overview of to which extent the pulses are varying over 
the pulse phase, we have calculated two values for each bin in the 
pulse profile: the modulation index and the R parameter. We define 
the modulation index as rrii = Oi/ fii, where Oi is the mean inten- 
sity in bin i and \ii is the standard deviation in the ith bin after 
integrating over the entire observation. While the modulation in- 
dex is a good indicator of whether there are persistent oscillations 
within the pulsed emission (such as e.g. drifting sub-pulses), it is 
not very sensitive to infrequent bursts of emission. In addition, it is 
hard to measure in observations with a low signal to noise value. To 
account also for non-persistent signal variation, we have calculated 
the R-parameter, R t = (MAX; — Hi)/oi, (where MAXi is the 
maximum intensity in the ith bin) as is described in Johnston et al. 
(2001). Ri indicates the presence of very bright, infrequent spikes 
of emission in the separate pulse bins. The off-pulse value for the 
R-parameter will increase with the number of rotations of the pul- 
sar due to Gaussian noise statistics, while the modulation index is 
undefined in the off-pulse region. 

Fig. 8 shows how the modulation index and the R-parameter is 
varying over the pulse profile for the 3. 1 GHz single pulse observa- 
tion. The minimum and maximum values of the modulation index 
for bins which are clearly within the on-pulse region are m m i„ = 
1.7 and m max = 4.5. These values are in the upper range of val- 
ues for ordinary pulsars (Weltevrede et al. 2005; Burke-Spolaor et 
al. submitted). From the R-parameter we can clearly see that the 
burst modulation is greater in the leading and the trailing edge of 
the on-pulse region than it is in the middle of the pulse. Over- 
all the R-parameter is high for this source, as well being among 
the upper range of measured values for a larger pulsar sample 
(Burke-Spolaor et al. submitted). This further demonstrates the ir- 
regular spikiness of the emission at each rotation. Also in the case 
of XTEJ1810-197 the integrated pulse profile consists of strong 
spiky sub-pulses (Serylak et al. 2009). However, some of the sub- 
pulses for XTEJ1810-197 could be considered giant pulses, but 
with broader pulse widths. Serylak et al. (2009) also report on a 
modulation index for the single pulses that is high on average and 
that increases with increasing observing frequency but varies be- 
tween components also within the same observing frequency. Sim- 
ilar spiky emission has been observed also in ordinary pulsars, 



such as B0656+14 (Weltevrede et al. 2006). Again in contrast to 
PSR J1622^1950, this source also emits giant pulses and here the 
stronger pulses appear spikier than the weaker ones. 



7 CONCLUSIONS 

The peak flux density of PSR J 1 622-4950 has decreased by a factor 
of ~2 since the discovery, and analyses of non-pulsed observations 
confirm this trend (Anderson et al. 2011b). The timing analysis of 
the source implies large variations in the rate of spin-down. We 
find a correlation between flux density decrease and spin-down de- 
crease, but are cautious not to over-interpret this result. 

The polarization is changing greatly between observations, but 
in general we observe a high degree of linear polarization and low 
circular polarization at frequencies higher than 1.4 GHz. The ob- 
servations at 1 .4 GHz are highly affected by interstellar scattering, 
which causes depolarization in the linear component, but not all 
changes in the linear polarization at this frequency can be due to 
propagation effects. We also see orthogonal phase jumps, flat PAs 
and changed handedness in the circular polarization in some of the 
observations at 1.4 GHz and also at higher observing frequencies. 
RVM fits imply that the geometry of PSR J1622-4950 is nearly 
aligned, with the line of sight remaining within the emission beam 
for large parts of the rotation. If this proves to be a general prop- 
erty for magnetars, it could provide insight in to why so few of the 
magnetars have observed radio pulsations. 

The single pulses from the magnetar are very narrow in com- 
parison to the width of the total integrated profile. The emission 
from each rotation consists of a few narrow spikes, that vary in 
longitude and separation between pulses. The widths of the sin- 
gle emission spikes appear to scale inversely with observed fre- 
quency, by getting narrower as the frequency gets higher. This 
could demonstrate that the emission spikes observed at different 
frequencies are emitted at different altitudes above the polar cap in 
the neutron star. The pulse flux distribution points towards a log- 
normal flux distribution of the spikes, without any signs of giant 
pulses from the magnetar. The spikiness in the emission results in 
a very high R-parameter across the pulse profile and a modulation 
index that is slowly increasing with pulse longitude and is higher 
in the pulse edges. 
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In all, PSR J1622-4950 has many properties that are very sim- 
ilar to the other radio magnetars, but also some things that differ, 
such as variations in linear polarization and PA swing on short time 
scales and a non-smooth frequency derivative evolution. We stress 
the importance of regular monitoring of these sources as the known 
sample is very small and the only chance to understand their emis- 
sion and their connection to other neutron stars is by continuing to 
observe their various properties. In particular, it is important to es- 
tablish whether the overall flux density is the only property that is 
changing on longer time scales, or if also other emission properties 
are varying. In order to discover more radio magnetars we need to 
understand the time scales of their on and off periods, and depend- 
ing on their special features, we may need to re-evalute the way 
searches for these sources should be performed. 
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